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Abstract 

The ionization of atomic hydrogen in intense laser fields is studied theoretically. The calculations 
were performed applying both quantummechanical and classical approaches. Treating the problem 
quantummechanically, the time dependent Schrodinger equation (TDSE) of our system was first 
transformed into a pseudo-momentum space and solved in this space iteratively. While neglecting 
the Coulomb potential during the solution of the TDSE we got the results in the Volkov approxima- 
tion, in the first order solution we taken into account the Coulomb potential as perturbation. The 
classical calculations were performed within the framework of the classical trajectory Monte-Carlo 
(CTMC) method. 

The double differential ionization probabilities are calculated for different laser pulses and a 
reasonable agreement was found between the theories. Major differences can be observed in the 
angular distribution of electrons at low electron energies between classical and the quantumme- 
chanical approaches. At high electron energies the differences disappear, which indicates that the 
generation of low energy electrons is of quantum type, and it is strongly influenced by the Coulomb 
potential, while the production of high energy electrons is of classical type and it is less influenced 
by the Coulomb interaction. Our results are also compared with the Coulomb- Volkov (CV) model 
calculations. 

PACS numbers: 32.80Fb, 42.50Hz 
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I. INTRODUCTION 



In the last years advanced laser facilities have achieved intensities of the order of 
lO^^W/cm'^ and pulse lengths of the order of 10 fs., which corresponds to few cycles of 
an electrical field of 800 nm wavelength l|, 121]. The process of interaction of such short and 
strong pulses with matter is a topic which has attracted much interest recently. 

II the multiphoton regime, many experimental d theoretical studies 



0,3 



11| have been performed. In the tunneling regime, on the other hand, recent 



experiments 12|, |l3| with linearly polarized lasers have shown structures in the momentum 
distribution of the photoionized electrons in rare gases which have not been fully understood 
so far. The role of interference in few-cycle pulses is actually been both theoretically and 



experimentally investigated 
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Moreover, the combined laser and charged particle induced electronic processes either 
in ion-atom or ion-surface collisions are of considerable interest both for basic and applied 



sciences p/7|. From the fundamental point of view they might broaden our general un- 
derstanding of the dynamics of atomic processes for laser-matter interactions and field-free 
collisions. As for the applications these studies can help us to find the way for the control of 
ultrashort quantum processes which are important in a number of applications, like in laser- 
driven fusion, in plasma heating, or in the development of fast optical electronic devices. The 
dynamics of atomic processes for the above mentioned interactions are not fully understood 
due to the lack of the exact and efficient theoretical models. For a detailed understanding 
of processes involved in the laser-matter interaction one needs to solve the time dependent 
Schrodinger equation (TDSE) for an atomic system in the radiation field, but its exact an- 



alytical solution is not known. Several 
these kinds of systems are known 18 



numerical solutions of the Schrodinger equation for 
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22l |. but they are very time consuming 



for large systems and at high radiation intensities they converge slowly. To overcome this 
problem there are several theoretical approaches, which are based on the simplification of 
the TDSE using different approximations depending on the laser field parameters. 

At low and moderate laser field intensities the time dependent perturbation theory 
(TDPT) is a well suited approximation for single- and multi-photon processes and even 
for the above threshold ionization [8]. In this case the TDSE is solved by considering the 
interaction between the laser field and the studied atomic system as a perturbation. The 
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TDPT approach breaks down at higher laser intensities when it fails to describe the "peak 
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suppression" in the above threshold ionization spectra 

At higher laser field intensities other non-perturbative processes emerge like high har- 
monics generation, tunneling ionization (TI), over the barrier ionization (OBI), which can 
not be described using TDPT and other approaches are necessary. 



Time independent, stationary approaches are developed based on the Floquet theory [25|. 
In the framework of the generalized Floquet theory a stationary treatment of a periodic 
or quasiperiodic TDSE is possible. The TDSE is reduced to a set of time-independent 
coupled equations or into a Floquet matrix eigenvalue problem. Approaches involving the 
Floquet theory are applied successfully to describe multiphoton excitation, above threshold 
ionization and dissociation, but their main disadvantage is that in most applications the 
analytical form of the involved Floquet states is not known and their numerical calculation 
is also time consuming. 

The n»st frequently used approaches are based ou the Ke.dysh theory Q. Q. The 
Keldysh theory is based on the assumption that on the final state wavefunction only the 
external laser field has a dominant infiuence and it can be considered as a momentum 
eigenstate. The main shortcoming of the Keldys theory is that it completely neglects the 
long ranged Coulomb interaction between the ionized electron and remaining target ion. 
The Coulomb interaction leads to phenomena like sub-peaks in the above threshold ioniza- 
tion spectra 28[, asymmetry in the spatial distribution of the ejected electron even for a 
symmetric few-cycle laser pulses 29|, which can not be explained in the simplified Keldysh 
formalism. There are two possible ways of including the Coulomb interaction in the Keldysh 
formalism. The first one is by making corrections in the transition matrix, and the second 
one is by making corrections in the Volkov wavefunction. Approaches based on the Keldysh 
theory using these corrections are applied with considerable success to study multiphoton 
and tunneling ionization of atomic systems 27|, 1301] . 

In the recent years the Coulomb- Volkov (CV) wavefunctions were used to describe pro- 



cesses in the presence of intense ultrashort laser field 31 



wavefunctions were introduced by Jain and Tzoar in 1978 
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331]. The Coulomb- Volkov 



3J to describe laser assisted col- 



lisions. Later on they were successfully applied to describe multiphoton and above threshold 



ionization 
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351 ]. The Coulomb- Volkov wavefunctions are also used to study ionization in 



tunneling and in over the barrier regime in the framework of sudden approximation 3l|, l32 ]. 
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This model is called Coulomb- Volkov model 

The accuracy of the results provided by the CV model is limited by the sudden approxi- 
mation (it provides accurate results only if the pulse duration is less then two orbital period 
of the active electron 32|) and by the fact that the CV wavefunctions are only an approx- 
imate solution of the TDSE for charged particle in the presence of external radiation field. 
From these two main limitations the first one is more restrictive, because it is limiting to 
the pulse duration in the attosecond region. Our goal is to construct a theoretical model, 
which provides results as good as or better than the CV model, but which does not have 
the limitation in pulse duration. 

Our approach is based on the approximate solution of the TDSE for quantum systems 
with one active electron, where the Coulomb interaction between the electron and the re- 
maining target ion can be considered as a perturbation during the external laser pulse. The 
time dependent wavefunction of the active electron is expanded in terms of Volkov wave- 
functions. The equation for the expansion coefficients obtained from the TDSE is solved 
with this approximation. 

The ionization process in the over the barrier regime is considered to be a classical one and 
it is believed that it can be described very well by classical models like classical trajectory 



Monte Carlo (CTMC) |20|, l32|, l36|. 



We present double differential ionization cross sections for different laser pulses using the 
present, CTMC and Volkov models. The ionization probability densities are presented as a 
function of the electron energy and ejection angle. 



In several works 37|, l38|, l39|] the effect of the Coulomb interaction on the ionization of 



the hydrogen atom by ultrashort laser pulses at different field intensities using different 
approaches is studied. In these studies the effect of the Coulomb potential is studied by 
employing different types of final state wavefunctions like the Coulomb- Volkov wavefunction 



39| , or first and second order Coulomb corrected Volkov wave 



'unctions (38|] , or by employing 
Coulomb corrections on the Keldysh- Faisal- Reiss theory 37(]. Moreover, the effect of the 



Coulom potential on the overall ionization process is also investigated. 

In the present paper we study the effect of the Coulomb potential during and after the 
laser pulse by analyzing the angular distribution of the electrons at given energies (extracted 
from the double differential cross section) and the ionization probability density (calculated 
from the double differential cross section by integration over the ejection angles) using various 
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model calculations. The dependence of the angular distribution of photoelectrons on the 
ejection energy allow us to investigate the transition between a quantum type and a classical 



type ionization also studied by Dimitrovski and Solov'ev 40| in the framework of the first 



Magnus approximation. Atomic units are used throughout the calculations. 

II. THEORY 

A. Characterization of the model 

The time evolution of atomic systems in the presence of one intense ultrashort laser pulse 
is investigated. The laser pulse, an external electromagnetic field from the point of view of 
the studied system, is defined by its electric component 

[ iEo siniut + 0o) sin2(^) if t G [0, r] 
I elsewhere 

where i is the polarization vector along which the laser pulse is linearly polarized, ut is the 
frequency of the carrier wave, 0o is the carrier envelope phase and r is the pulse duration. 
The carrier envelope phase is set as follows: 

leading to a time symmetric laser pulse. The laser pulse is composed of a carrier wave 
sin(u;t — ^ — f ), which is modulated by a sine-square envelope function sin^(^). Figure [T] 
shows the shape of the laser pulses used in our calculations. 

In the time evolution of the studied system one may distinguish three main time intervals. 

In the first time interval {t < 0) the laser field is not yet switched on and the studied 
system is in a field free eigenstate. The Hamilton operator of the system in these conditions 
is 

Hi = ^ + V{r-), (3) 

where ^(r) is the Coulomb potential between the active electron and the rest of the system. 
The eigenvectors and eigenfunctions of this Hamilton operator are considered to be known 

Hjij, = E.ij, (4) 
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FIG. 1: (Color online) The electric components of the laser pulses used in our calculations. Solid 
line - r = 3 a.u.; Dotted line - r = 5 a.u.; Dashed line - r = 10 a.u. 

and the initial state of our system coincides with eigenstate V'j- Under these conditions the 
time dependent wavefunction of the system is 

vl>,(t) = V'.e-^^'*. (5) 

During the second time interval (0 < t < r) the laser field is switched on and the Hamilton 
operator of the system becomes 

Hn = Y + rE + V{r), (6) 

where fE is the interaction term between the laser field and the active electron expressed 
in E gauge (Electric field gauge, length gauge). The dipole approximation was implicitly 
applied when the spatial dependence of E was neglected in Eq. ([T]). The time dependent 
Schrodinger equation 

t^^^nit)=Hn^nit), (7) 

with \l///(t) time dependent wavefunction of the active electron in the presence of the laser 
pulse, does not have known analytical solution. In order to determine the time evolution of 
the studied system one needs to know the \l//7(t) wavefunction. In the present approach the 
wavefunction is considered in the following form 

■^n{t) = j dkc{k,t)-<^v{k,t), (8) 



where are the Volkov wavef unctions. The Volkov wavefunctions are the solutions 

of the TDSE in dipole approximation for a charged particle in radiation field 

ii-^vik,t) = (^ + fE] ^vik,t), (9) 



dt ' ' \2 
and they can be expressed as 

In the above expression 

A{t) = - I E{t')dt' (11) 
Jo 

is the vector potential of the electromagnetic field. 

The time dependent wavefunction is well defined by the expansion coefficients c{k, t), so if 
one is interested in the time evolution of the system is enough to calculate these coefficients. 

By substituting the time dependent wavefunction of Eq. into the TDSE given by Eq. 
(I?]) can be obtained 

J dkc{k,t)^vik,t) = (^^ + V{r)+fE^ J dkc{k,t)mv{k,t). (12) 
Using Eq. ([9]) it is easy to prove that 

j dkc{k,t)i-^^vik,t) = (^y + ^-^^ J dkc{k,t)^v{k,t), (13) 
which can be used to simplify Eq. ffT^ 

i J dk^v{k,t)^^c{k,t) = j dkc{k,t)V{r)^vik,t). (14) 



^ j dk( ^c{k,t) ] e~UJ^*'(^>^>'))^ / dfe'Ck^P+'^mr 



Eq. fll4p can be converted into a more favorable form by transforming it into momentum 
space 

i I dh I , 

j dkc{k, t)e-i 'it'ik+Anr J rfr-y (r)e*(^^-^^+^>))n (15) 
After basic mathematical operations and the substitution q = p—A{t), this equation becomes 

,(2vr)3^c(g,t)e-i/o'^*'(^+^V))^ = f rf^c(^,t)e-^/o'i*'(MO)^ f drV{f)e'^'^-P+^^'^>^. (16) 
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By introducing the notation s = k — q and by restructuring Eq. ffTUj) for the expansion 
coefficients c{k, t) can be obtained the following integro-differential equation 

t) = -^^l dsc{s + q, t)e-yo <^*'< W2A(0) J drV{r)e''\ (17) 

We note, however, that Eq. ( ITTl) is equivalent with the Schrodinger equation given by 
Eq. ([7]). By solving Eq. (1171) one obtains directly the time dependent wavefunction in 
momentum space, which carries all the information about the studied system. Eq. ffT7|) can 
be solved numerically, but this direct approach needs large computational resources, and in 
several cases it is not worth to do it, because by introducing some approximations one may 
keep almost the same accuracy. 

In the third time interval (t > r) the laser field is switched off and the Hamilton operator 
of the system can be expressed as 

Hin = ^ + V{r), (18) 
which is identical with the one given in the first time interval. 

The time dependent wavefunction of the system in this time interval can be given as a linear 
combination of stationary-state wavefunctions 

^m(r, t) = J2 ^^e-'^"* + [ dkcf{k),l:fe-'^f\ (19) 

b 

where V'fe represents bound, while represents free states. For convenience the free states 
are represented by plane waves, so we have 

iPf = e'^'' (20) 

and 

k"^ 

Ef = y. (21) 

One of the basic properties of the time dependent wavefunction, which describes the 
evolution of a real system is that it is continuous over time. 
From the continuity condition at time t = one obtains 

^i{t = 0) = ^n{t = 0), (22) 

which has the following explicit form 

= j dkc{k,t = 0)e'^''. (23) 
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From this expression one may get the initial condition for Eq. (1171) 

c(g,t = 0) = ^(e^^"^|V'.) (24) 

It is worth to notice that the above initial condition is in fact the initial state wavefunction 
in the momentum space, because at t = the Volkov wavefunctions are reduced to plane 
waves. 

From continuity condition at t = r we can write: 

^n{t = r) = ^niit = r), (25) 

which can be given explicitly as 

f dkcik, T)e-y^ dt'Ck+At'))\^riUA(r)) ^ ^ ^^^-iE,r ^ f dkCf{k)e''~^'^-'^f^\ (26) 
J b ^ 

We are interested in the transition probability to a free final state, which means that 
we need to know Cf{k). This can be obtained by transforming the above equation into 
momentum space and by neglecting the contribution of bond states 

Cf{p) = c{p- l(r), ^)e-^[^(p-^(-).-)-^/-l, (27) 
where we used the notation 

E{k,T) = \j^ dt\k + A{t')f. (28) 

The Cf{'p) given by Eq. ([27]) is the final-state wavefunction in momentum space. 

The transition probability from the initial state ipi to a free final state ipf, with a well defined 

momentum p, is given as 

P,^f = {27rr\cf{p)f = i2nf c(p-l(r),r) (29) 



B. Volkov model as a zero order approximation 

The critical point of our model is how to solve Eq. ( ITTl) or in other words how to 
propagate our system in the presence of external laser field in the second time interval. 
The simplest way is by neglecting completely the Coulomb potential {V{f) = 0). This 
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approximation provides good results only for high laser field intensities. In this framework 
Eq. ffTTl) becomes 

^c(g,t) = 0, (30) 
which has the following analytical solution using the initial condition (Elj) 

c(g,t)^cW(g-) = ^(e^nV^.). (31) 

From Eq. (|271) the final state wavefunction in momentum space can be written as: 

Cf{p) = c\p- l(r))e-^[^(P-^W'^)-^/^, (32) 

which is the initial state wavefunction in the momentum space shifted by the momentum 
transfer A{t) gained from the external laser field. 

C. First order approximation 

In most cases the Volkov model (see Eq. ( !32l) ) does not provide accurate results, because 
the Coulomb interaction at moderate intensities can not be totally neglected. At moderate 
laser intensities one can assume that the infiuence of the Coulomb interaction on the evolu- 
tion of the system in the second time interval is small and the expansion coefficients c{k^ t) 
are close to that ones provided by the Volkov model flH2|) . Based on this arguments Eq. ffTT]) 
can be simplified by replacing c{k,t) in the right hand side by c^^\k) as follows 

|-cW(g,t) = [ dscW(s + g-)e-^^o'^*'<W2AV)) f ^^v{f)e^^. (33) 

dt (27r)^ J J 

The same equation can be obtained using a different approach by considering the Coulomb 
interaction as a small perturbation and by keeping only the first order terms in V{r). The 
advantage of this approximation is that it eliminates the direct coupling between the expan- 
sion coefficients c{k,t), which makes easier and faster the solution of Eq. ( |T71) . 
Eq. (l33i) can be simplified and its solution can be given as simple integral 

c«(g, t) = cW(g-) -J^,[ dt'Iiq, t'), (34) 



where 



I{q,t) = j rfsc(0)(s + g-)e-no'^*'<W2A(t')) j ^^y^^^^sr_ (35) 
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D. The Classical Trajectory Monte Carlo simulation 



In the last two decades there was a great revival of the classical trajectory Monte Car 



(CTMC) calculations applied in atomic collisions involving three or more particles |41 |. 
This approximation seems to be useful in treating atomic collisions where the quantum 
mechanical ones become very complicated or unfeasible. This is the case usually when 
higher order perturbations should be applied or many particles take part in the processes. 
The CTMC method has been quite successful in dealing with the ionization process in laser- 
atom collisions, when, instead of the charged particles, electromagnetic fields are used for 
excitation of the target. The CTMC method is a nonperturbativ method, where classical 
equations of motions are solved numerically. A micro-canonical ensemble characterizes the 
initial state of the target. In this work, the initial conditions of the target are taken from 
this ensemble, which is constrained to an initial binding energy of H(ls) (0.5 a.u.). 

In the present CTMC approach, Newton's classical nonrelativistic equations of motions 
are solved 



42 



43 



44l | numerically when an external laser field given by Eq. ([T]) is included. 
For the given initial parameters, Newton's equations of motion were integrated with respect 
to time as an independent variable by the standard Runge-Kutta method until the real 
exit channels are obtained. For the ionization channel the final energy and the scattering 
angles (polar and azimuth) of the projectile and the ionized electron were recorded. These 
parameters were calculated at large separation of the ionized electron and the target nucleus, 
where the Coulomb interaction is negligible. 

The single and double differential ionization probabilities {Pi) were computed with the 
following formulas: 

(36) 



dE NAE 

dPi N,, 



dn NAn ■ 

d'P 



dEdn NAEAQ 
The standard deviation for a differential probabilities is defined through: 



(37) 



(38) 
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AP, = P, ^ . (39) 

[N-Ni\ ^ ' 

In Eqs. (1361 - is the total number of trajectories calculated for the given colli- 

sion system, Ni is the number of trajectories that satisfy the criteria for ionization under 
consideration in the energy interval AE and the emission angle interval AQ of the electron. 

III. IONIZATION OF THE HYDROGEN ATOM 

The theoretical approaches presented above are applied to describe the ionization of the 
hydrogen atom in over the barrier regime. We choose this system because the calculations 
are easy to perform and because there are several theoretical studies on this system 
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39j and it is easy to compare our results with others found in the literature |32| 
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Using the Is orbital of the hydrogen atom as initial state wavefunction 



ipi = ^e-' (40) 



one obtains the following initial condition for Eq. ([T 



= , ^^,^-^^2 ■ (41) 

Using Eq. (HTi) the ionization probability in Volkov approximation can be expressed as 

Pi->f{p) = ■ (42) 

vr(l + (p-A(r))2)4 

The first order model can also be adapted very easily by using the Coulomb potential 

V{r) = -^ (43) 
in Eq.( l35ll . which can be significantly simplified by performing some of the involved integrals 



^ ~ r_ ™+c» Jo Jo (1 + s2 + g2 + 2sq cos(^,))2 ^ 
^_^(,2i+2s,cos(0.)t+2/Wscos(e.)cos(e.)) j^(2s/(t) sin(^,) sin(e,))(cos(sw) - 1), (44) 

where Jo is the Bessel function of the first kind. The angle between s and g is Og, while the 
angle between e and q is 9^ and 

fit) = [ A{t')dt'. (45) 



By substituting Eq. (H^ into Eq. (!34l) and by performing the remaining integrals numeri- 
cally one obtains the expansion coefficient in first order approximation. From the obtained 
expansion coefficients the ionization probability is calculated by using Eq.( 
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IV. RESULTS AND DISCUSSION 



Calculations are performed using laser pulses with duration r of 3 a.u., 5 a.u. and 10 a.u. 
at two different field intensities {Eq = 1 a.u. and Eq = 10 a.u) for each pulse duration. The 
shape of the pulses can be seen on FIG. [H 

The double differential ionization cross sections calculated using the Volkov, the first order 
and the CTMC model are presented on Figure [21 where the ionization probability density 
is plotted as a function of the electron energy and the ejection angle. The ejection angle 
is measured from the polarization vector i, which is in the direction of z axis. Due to the 
spherical symmetry of the ground state wavefunction and of the Coulomb potential the 
double differential ionization cross section has a rotational symmetry, where the rotational 
axis is in the direction of the polarization vector e of the laser field. Due to this symmetry 
the yOz positive semi-plane carries all the information about the ionization cross section. 
At first sight one can observe that at a large scale all three models predicts the same cross 
sections. In each approach the electrons are ejected with maximum probability along the 
polarization vector i with energy around the value , which is gained by the momentum 
transfer A{t) from the external laser field. 

After detailed analysis, important differences can be observed. In the case of the first 
order and CTMC model the maxima of the predicted cross sections are shifted toward smaller 
energies. This shift is caused by the Coulomb attraction during the ionization between the 
active electron and the rest of the system, in a classical picture the ejected electrons are 
decelerated by the Coulomb attraction. 

Other important differences can be identified at the angular distribution of the ejected 
electrons at fixed energies (We) as shown in Figure El The maximum value of each angular 
distribution is normalized to unity to allow an easier comparison. 

At the low part of the energy spectrum (see FIG. El (a) and (d), with below the 
peak observable on FIG. 2) significant differences are observed between the predicted first 
order and CTMC angular distributions. The observed differences implies the quantum 
nature of the ionization for low energy electrons, which can not be described correctly by 
classical calculations. These differences start to disappear at higher electron energies (see 
FIG. O (b) and (e)), where both distributions roughly are the same. However some minor 
differences still exist showing a transition between the quantum and classical nature of the 
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FIG. 2: (Color online) Two-dimensional ionization probability density in the yOz positive semi- 
plane as a function of the electron energy and the ejection angle. First column - Volkov results; 
Second column - first order results; Third column - CTMC results; First row - r = 3 a.u. and 
£"0 = 1 a.u.; Second row - r = 3 a.u. and £"0 = 10 a.u.; Third row - r = 5 a.u. and Eq = 1 a.u.; 
Fourth row - r = 10 a.u. and Eq = 1 a.u. 
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FIG. 3: (Color online) Angular distribution of the photoelectrons for given electron energies (We). 
Solid line - first order; Dotted line - Volkov; Squares - CTMC; First column - r = 5 a.u. and 
£^0 = 1 a-u.; Second column - r = 10 a.u. and Eq = 1 a.u.; (a) We = 0.02 a.u.; (b) We = 0.125 
a.u.; (c) We = 3.5 a.u.; (d) We = 4.5 a.u.; (e) We = 6 a.u.; (f) We = 10 a.u. 
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ionization. At high energies, the differences observed at lower energies completely disappear 
(see FIG. E] (e) and (f )) and the distributions predicted by first order and CTMC model are 
in good agreement indicating that the ionization mechanism for electrons ejected at high 
energies behave classically. This evolution of the cross section and the transition between the 
quantum and classical type ionization was also studied by Dimitrovski and Solov'ev 401] in 
the framework of first Magnus approximation. They obtained the same transition between 
quantum and classical type ionization by varying the net momentum transfer. 

At the low part of the energy spectrum (see FIG. [H] (a) and (d)) large discrepancies are 
also observed between the predicted first order and Volkov angular distributions. These 
discrepancies show the significant influence of the Coulomb interaction on the ionization of 
the low energy electrons. This disagreement between the discussed models diminishes and 
completely disappears at higher electron energies (see FIG. E] (b, c, e, f)), which indicates 
that the Coulomb interaction has less influence on the ionization of the high energy electrons. 
FigureHlshows the net angular distributions of the ejected electrons obtained from the double 
differential cross section by integrating over the electron energies for different pulse durations. 
The angular distributions predicted by different models are in good agreement with each 
other, because the contribution of the low energy electrons in the net angular distribution 
is negligible, so the influence of the Coulomb potential and the quantum nature of the 
ionization can be also neglected. On Figure Owe have represented the ionization probability 
densities obtained from the double differential ionization cross sections by integrating over 
the ejection angle. The probability densities are calculated using the first order, the Volkov 
and CTMC models. Figure [5] also shows the TDSE and CV results calculated by Duchateau 
et. al. 32|. At high laser field intensity {Eq = 10 a.u.) a very good agreement is observed 



between TDSE and CTMC results on Figure (b). At lower laser intensity {Eq = 1 a.u.) 
the agreement is acceptable (see FIG.E](b) and (c) ), but not as good as in the previous case, 
because at these intensities tunneling ionization also takes place, which is not considered in 
our CTMC model. This good agreement between classical and quantum approaches at high 
laser field intensities was also confirmed by several other studies 20|, l32| , which proves the 
classical nature of the over the barrier ionization. 

The main difference between the Volkov and the first order model is that the Volkov 
model completely neglects the Coulomb attraction during the ionization process. The effect 
of the Coulomb interaction on the photoelectrons during the ionization can be studied by 
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FIG. 4: (Color online) The net angular distribution of the photoelectrons. Solid line - first order; 
Dotted line - Volkov; Squares - CTMC; (a) - r = 5 a.u. and Eq = 1 a.u.; (b) - r = 10 a.u. and 
£"0 = 1 a.u. 

comparing the Volkov and first order ionization probability densities. At the high energy part 
of each photoelectron spectrum presented on Figure El a good agreement can be observed 
between the Volkov and first order model. This agreement breaks down at the lower energy 
part of the spectra, where the first order spectrum is shifted toward the lower energies. The 
above behavior was also observed in the angular distribution of ionized electrons and it can 
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FIG. 5: (Color online) Ionization probability density as function of the electron energy. Solid line 



32]; Dash-dotted line 



- first order; Dotted line - Volkov; Squares - CTMC; Dashed line - TDSE 

n 

CV [32] First column - Eq = 1 a.u.; Second column - Eq = 10 a.u.; First row - r = 3 a.u.; Second 
row - T = 5 a.u.; Third row - r = 10 a.u. 
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FIG. 6: (Color online) Two possible electron trajectories during the ionization corresponding to 
photoelectrons with high (Solid line) and low (Dotted line) ejection energies. The black dot is the 
residual proton and the net momentum transfer is along the z axis. 

be explained using a very simple intuitive picture presented on Figure [HI The electrons 
with high ejection energy have their momentum in the initial state in the same direction 
as the net momentum transfer. In this way the trajectory of the high energy electrons 
leads directly away from the core, with a very small portion close to the core where the 
Coulomb potential has a significant influence. On the other hand the electrons with lower 
ejection energies have a momentum in the initial state, which leads in the opposite direction 
of the net impulse transfer. In this way the low energy electrons need to "go around" the 
core, leading to a trajectory, which has a long portion close to the core, where they can be 
influenced significantly by the Coulomb interaction. 



Results obtained by Duchateau et. al. 



32| using Coulomb- Volkov wavefunctions and the 



sudden approximation is also included in our analysis. In the CV model the influence of 
the Coulomb interaction on the ionization process is not included, it is taken into account 
only in the final state wavefunction. The comparison of CTMC, CV and first order results 
allow us to study in detail the influence of the Coulomb interaction during and after the 
laser pulse. 

At low intensities for short pulse, the agreement between the CV and CTMC results 
is better than between the CTMC and first order results (see FIG. (a)). For this pulse 
parameters the net momentum transfer is small leading to photoelectrons with small mo- 



19 



mentum, which can not be described correctly by the plane waves employed by the first 
order model. 

At the same intensity but for longer pulses, major discrepancies can be observed between 
CV and CTMC results (see FIG. |5] (c, e)), because for these pulse durations the sudden 
approximation no longer applies. The same discrepancies are observed between first order 
and CTMC results, because the final state wavefunction is still inaccurate. 

At high laser field intensities the agreement between the first order and CTMC results is 
better than the agreement between CV and CTMC results (see FIG. El (b)). In this case the 
momentum transfer is high leading to photoelectrons with high energies, where the plane 
wave approximation for the final state is accurate. At the low energy part of the spectra 
presented on Figure [5] (b, d, f) the first order results show good agreement with the CTMC 
results (considered identical with the TDSE results at this high intensity), but discrepancies 
exists at the high energy part of the spectra. Both first order and CV models predict a higher 
ionization probability than the CTMC model at the high energy part of the spectra, which 
in the case of the first order model can be corrected by taking into account the Coulomb 
interaction after the laser pulse. 



V. CONCLUSIONS AND OUTLOOK 



In the present work a general approach for the ionization of atomic systems by ultra- 
short laser pulses in one active electron approximation was presented. The time dependent 
Schrodinger equation of the studied system was transformed into a pseudo-momentum space, 
where it was solved in an iterative way. In the zero order solution the Coulomb potential 
was neglected (Volkov model), while in the first order solution it was taken into account 
during the ionization as a perturbation. 

Calculations were performed for the ionization of the hydrogen atom in over the barrier 
regime using the Volkov and the first order solution of the TDSE. Classical trajectory Monte- 
Carlo calculations were also performed. The obtained results were analyzed and compared 



3^ using TDSE and CV model. 



with results obtained by Duchateau and coworkers 

The double differential ionization cross sections were calculated using Volkov, first order 
and CTMC models for different laser pulses. A good agreement was found between the 
results using these three models at high laser field intensities. This good agreement between 
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classical and quantum calculations was also confirmed by other groups [20|, l32|. At lower 
intensities, however, small discrepancies appeared due to the tunneling ionization. 

The first order and CTMC results were shifted toward smaller energies due to the 
Coulomb attraction of the remaining target ion. 

More differences were identified by analyzing the angular distribution of the ejected elec- 
trons. For low energy electrons major discrepancies have been found between CTMC and 
first order results, because the ionization of this low energy electrons is a quantum process, 
which can not be described in a classical approach. At higher electron energies, the observed 
discrepancies disappear leading to the conclusion that the ionization process of the high en- 
ergy electrons can be described classically. At low energies, differences were also observed 
between the first order and Volkov results, which disappeared at higher electron energies 
indicating that during the ionization electrons with higher energies are less influenced by the 
Coulomb interaction. The same observation was made by analyzing the ionization proba- 
bility densities, where significant differences between the Volkov and first order results were 
observed for low energy photoelectrons. A simple explanation for this behavior was found by 
considering the possible classical trajectories of the electrons during the ionization. Similar 



conclusion was also drawn by Zhang et. al. 



39j in their study, where they found that the 



influence of the Coulomb interaction depends mainly on the energy of the photoelectrons 
rather than on the laser field intensities. 

The influence of the Coulomb potential during and after the laser pulse was studied. For 
short laser pulses with low momentum transfer we found that the Coulomb potential has 
an important role after the laser pulse is switched off. In the case of longer pulses, with 
low momentum transfer, the Coulomb potential is equally important during and after the 
pulse, while in the case of high momentum transfer the Coulomb interaction has a greater 
influence during the laser pulse. 

A good agreement between first order and CTMC results was found at high laser field 
intensities, where the momentum was high, while at lower intensities with low momentum 
the agreement was acceptable, comparable with the agreement between the CV and CTMC 
results. 

Our model can be further improved by considering the influence of the Coulomb potential 
on the final state wavefunction, by using Coulomb wavefunctions instead of simple plane 
waves for the final state, which will be the subject of our future investigations. 
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